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Abstract 

We study the influence of turbulent viscosity variability on the proper- 
ties of self-propelled body wake model. In addition to the already known 
integrals of motion obtained with constant turbulent viscosity, we obtain 
new ones. The presence of new integrals of motion leads, in particular, to 
changes in the behavior of the width and profile of the wake leading to its 
conservation. 

Study the behavior of long-lived turbulent wakes behind moving bodies is not 
only important from a scientific [3J GH H3] but also practical point of view [H [71 
[T0| . On the one hand the evolution of such wakes can affect the operation of 
airports, and on the other - mixing processes in the surface layer of the sea or 
ocean. The turbulent wakes behind ships have attracted special interest, largely 
due to the need to interpret the radar observations of ship wakes [51 [TU1 HU H2] . 

One of the most important models of the turbulent wakes behavior, from a 
practical point of view, is the model of self-propelled body wake [U [51 El [Zl [H] > 
and the main characteristic of this wake is turbulent viscosity,^ . Analytical 
approaches so far assumed a constant viscosity v t — const , which leads to 
the wake width at the distance x equal to W(x) ~ a;" (where a — 1/5 for 
an axisymmetric wake and for a = 1/4 a plane one). However, in real wakes, 
turbulent viscosity may be a function of the coordinates. In this paper we study 
the consequence of this dependence. 

Consider the axisymmetric wake in which the mean flow in planes passing 
through the axis of symmetry is identical and weakly depends on the distance 
x. Simultaneously, we examine flows whose principal mean velocity component 
is directed along the x-axis. This means that radial component of the mean 
velocity \ U r \<<\ U x \. If L is a typical spatial scale of inhomogeneity in 
the direction of the wake (x -direction) while I is a typical spatial scale of 
inhomogeneity in the transverse direction (r), I « i(in this paper we use the 
notations in conformity with |13|). Then we denote by the mean flow velocity 
in x direction, and by U s the maximum value of | Uq — U x \ (U s << Uq far 
from the body, ) of the cross-wake variation of the mean velocity component in 
direction. Following [13J , an equation for the momentum transfer in x direction 
can be written as: 

Ox r Or 

where U — U x — Uo , q xr =< u x u r >, Ui is a turbulent component of the velocity 
field, and <...> denotes averaging over an ensemble. Taking into account self- 
preservation hypothesis,: the velocity defect and the Reynolds stress become 
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invariant with respect to x, and they are expressed in terms of the local length 
and velocity scales l(x) and U s {x) . That is. jj- — f(r/l) and = —g(r/l). 
We substitute these expressions into the equation of motion fTJ) : 

U l dU s U dl > , a 

-uf^ f+ ir s TJ f = 9 + ? (2) 

where £ = r/l(x). To satisfy this equation for all x we obtain / = Ax a , U s — 
Bx 01-1 , where A and B are constants. Substitution into ^ gives 

fl (a-l)/ -«*£/' = —(&) (3) 

where s = (Uq/B). The functions /(£) and g(£) are still to be found. 
Since the turbulent viscosity is defined by the relation 

-q X r = v t —, (4) 
the equation of motion acquires the following form: 

Uo ^L = I A (5) 

where i/ t = Bx 2a ^ 1 (g j '/'), and prime denotes £ derivative. 

oo 

The integral Iq = J Urdr is proportional to the total momentum across 
o 

oo oo 

the wake and is conserved. Indeed, 4^ J UoUrdr = — f -§p (rq X r) dr = 0. The 

o o 

conservation of momentum, together with the assumption that U oc x/(£), 
immediately gives Iq oc a; 3 " -1 J Q °° an d therefore we obtain a well-known 

value a = 1/3. 

However, it fails for self-propelled body wake, where the total momentum 
vanishes (Iq = 0) . In this case, multiplying by r m+1 and assuming the conver- 

oo 

gence, one has = - J Q °° r m J^(rq xr )dr , where I m = j U Ur m+1 dr. Using 

o 

turbulent viscosity, we obtain: 

*u - /> » Up, * 

ax 7 or \ or I 

o v ' 

Integrating in parts twice and taking into account that U — > and ^ — >• as 

r — > oo : 

oo 



Tx Im=m [u§;;(r m vt)dr 
o 



If -7^(r m v t ) oc r then J £/J^ (r m v t ) dr (x Iq = Oand 7 m is conserved, 
o 
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Figure 1: Wake profiles for certain parameters m. m = 2 corresponds to 
v t — const. 



In this case, v t = K(x)^ 2 m , where K(x) is some function on x. Conserva- 
tion of I m , together with self-preservation hypothesis, gives 

a =zA^ ( 6 ) 



It should be noted that similar results can be obtained for the plane flow. In 

this case we obtain a = — V^r . 

Thus, with turbulent viscosity definition and conserved I m , we obtain 
K(x)£ 2 - rn = Bx 2a - 1 g/f, and therefore K{x) = Bx 2a ~ x and g = £ 2 ~ m f' 

.The latter relation can be now substituted into ([3]). Eventually, we arrive at 

the following equation: 
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r + ^±^r- 2 / - o ( 7) 



The solution of this equation for certain values of m (and s = 1) is shown in 
Fig. 1. 

It is worth mentioning that although the turbulence viscosity coefficient, 
formally defined by ((H) , diverges as £ — ► , physical quantities remain well- de- 
fined. Really, U(0) is constant and nonzero, near £ = : / oc £ i m ^ 1 exp(— a£ m ) 
(where a is a constant) , and q xr oc g oc £, 2 ~ m f oc £ . 

As we can see in Fig. 1, the growth of m (more rapid decrease of the 
turbulent viscosity with £) leads to an effective reduction of the wake width. At 
the same time the zone characterized by U /U s — const increases. 

Thus, in the present work we have suggested to get rid of the assumption 
of the uniform turbulent viscosity for a self-propelled body wake. Instead, we 
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consider a power- law dependence vt oc r 2 ~ m , which allows us to construct a new 
integral of motion, I m = J °° r m Urdr instead of the vanishing I = J °° Urdr. 
Each of these integrals leads to the formation of a corresponding flow profile. 
Simultaneously, the increase in the parameter m leads to an increase in the flow 
"core" and slower growth of its width. 
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